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Abstract 

' The Hartree-Fock ground-state phase diagram of the one-dimensional 

, Hubbard model is calculated in the fi — U plane, restricted to phases with 

f"^ ■ no charge density modulation. This allows antiferromagnetism, saturated 

\^ ' ferromagnetism, spiral spin density waves and a collinear structure with 

■ unit cell Till- The spiral phase is unstable against phase separation near 

quarter-, half- and three-quarter-filling. For large U this occurs at hole 
(or electron) doping of (3t/7r 2 [/) 1//3 from half filling. 
PACS numbers 75.10Lp, 75.25+z 

^ \ 

^ ■ 1 Introduction 

Although the exact ground state of the one-dimensional one-band Hubbard 
model with nearest-neighbour hopping is well known|l|, systematic Hartree- 
, Fock (hf) studies of the same model are still of value. They are correct in the 

!_h ' limit of large degeneracy, they represent an effective Hamiltonian in functional 

integral formalism, and they provide a toy model for understanding the magnetic 
phases of two- and three-dimensional systems (such as cuprate superconductors 
and transition metal alloys). Such studies also raise suggestive connections with 
state-selection problems in frustrated Heisenberg magnets. 

The present author has obtained the hf ground state phase diagram of the 
one-dimensional nearest-neighbour Hubbard model [^) 
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for arbitrary band filling n (0 < n < 2). The calculation was restricted to phases 
without charge density modulation; subject to this restriction, the full n — U 
phase diagram was obtained. 
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The work here extends these results to the grand canonical case, presenting 
the ji—U phase diagram and the asymptotic form of the phase boundaries. Full 
details of the theory are given in the above reference. 

2 Hartree-Fock phase diagram 

The unrestricted hf approximation minimises the expectation value of the Hub- 
bard Hamiltonian ([j]) in the space of Slater determinants. These states are 
ground states of a non-interacting many-electron system in a potential specified 
by variational parameters. 

We restrict consideration to the uniform phases, where the only spatial de- 
pendence is in the local magnetisation directions, and to macroscopic phase 
separation. In the latter case two uniform phases are in equilibrium, separated 
by a domain wall. There are two families of uniform phases: 

• SSDW {spiral spin density wave, j /* — »\,) of continuously varying wave 
vector Q. This has limiting cases 

Q = 0: fm (saturated ferromagnetism, tttt)- 
Q = 7r: AFM (antiferromagnetism, Titi)- 

• dsdw (double spin density wave, or two interpenetrating antiferromagnctic 
sublattices with Neel vectors canted at an angle 9, 

For each point (n, U) the energy of each family is minimised with respect to 
the exchange splitting and the angle (Q or 9). The hf energy _Bhf(«, U) is the 
lower of the minima for ssdw and dsdw. 

The condition for stability against macroscopic phase separation is that 
E-iw(n, U) be a convex function of n, or equivalently that the chemical potential 



be an increasing function of n. If this does not hold, a Maxwell construction |3[, 
shown in figure 0(a), determines the fraction of each phase. 

The figure shows the generic picture for U > 8.7t. For n\ < n < 1 (and 
1 < n < 2 — m) the hf ground state consists of a uniform afm phase with n = 1 
in equilibrium with a hole-rich (electron-rich) uniform fm phase with n = nx 
(n = 2 — ni). The chemical potential is pinned within the afm gap and within 
the lower (upper) band of the fm phase. Additional carriers will simply move 
the domain wall. A similar effect is seen in supercell calculations of collinear 
configurations ||. For n < n\ (and 2 — n\ < n) the ground state is fm. Thus the 
ssdw is always unstable against separation into afm and fm phases. (Note that 
the ssdw energy peels off from the fm energy, indicating Q > 0, only above n\ .) 
The situation is slightly different for U < 8.7t: near n = 1 the ssdw is unstable 
against phase separation into an afm phase and a ssdw of longer wavelength. 

The previous paperQ presented the phase diagram in the n — U plane. 
Figure 0(b) shows the phase diagram in the fi — U plane. The areas shown 
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Figure 1: (a) The hf energies for U = lOt. Bold curve: ssdw energy; dot- 
ted curve: fm energy; thin line: energy of phase-separated state, (b) The hf 
phase diagram of the Hubbard model. The afm phase has filling n = 1, and 
the collinear dsdw phases (denoted 'D') have filling n = 1/2,3/2. The region 
marked 'full' corresponds to filled bands (n — 2). 



here are pure phases, separated by second-order transitions (between fm and 
finite-Q ssdw) and first-order transitions (between ssdw and collinear dsdw, 
between ssdw and afm, and between fm and afm). For large U, as a function 
of chemical potential there are two discontinuities in filling factor and in wave 
number. As U is decreased, further discontinuities appear (although only the 
dsdw phase has been calculated). 

3 Form of phase boundaries 

We can compute the asymptotic form of the fm-afm phase boundary for large 
U . By expanding the hf energy of the afm state at n = 1 we find the gap to be 

A = U - At 2 /U + At 4 /U 3 + ■■■ (3) 

and the energy 

E KF (1,U) = -2t 2 /U + 2t 4 /U 3 + ■■■. (4) 

(This of course exceeds the Bethe Ansatz ground state energy — 41n2t 2 /C/Q.) 
Drawing a tangent to the fm energy —(2t/ir) sinn7T gives the form of the phase 
boundary for large U as 

"' = 1 -(^) V3 -IoV + "- < 5 » 
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or 
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4 Discussion 



The full hf solution to the one-dimensional Hubbard model is more difficult 
to compute, in general requiring a 4A a -dimcnsional minimisation. One might 
speculate on a possible devil's staircase for small U: tongues, such as that 
for dsdw, would exist for all rational fillings. Similar behaviour is seen in the 
Falicov-Kimball model (where one spin state is immobile) || . For larger U (in the 
Hubbard model) the tongues would disappear, thinning out the staircase. We 
would indeed expect the uniform ssdw phases to distort in such a way as to open 
a gap at the Fermi surface; the dsdw can be seen in such a wayQ. The phase 
separation seen here may appear in a microscopic form (as a soliton lattice) 
rather than the macroscopic form discussed here. Realistic terms, absent in the 
one-band Hubbard model, that suppress long-wavelength charge fluctuations 
would tend to prevent macroscopic phase separation. 

It is also interesting to note connections with the problem of state selection 
in frustrated Heisenberg magnets, dsdw phases with varying 9 are strictly 
degenerate in a classical Heisenberg model; here it is the itinerant nature of the 
magnetism that selects the collinear state. 

We end by confessing that these results bear little relation to the true ground 
state of the one-dimensional one-band Hubbard model, which does not break 
symmetry in this way. However, they may provide a useful starting point for the 
Hubbard model in two and three dimensions and many-band Hubbard models 
in one dimension. 
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